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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2^36 

HEAT DELIVERY IN A COMPRESSIBLE FLOW AND 
APPLICATIONS TO HOT-WIRE ANEMOMETRY 


By Chan-Mou Tchen 


SUMMARY 


In a two- dimensional field a generalized potential theory applicable 
to nonadiabatic and rotational flow is developed. Three partial differ- 
ential equations are first obtained determining the three variables which 
are: Distribution of additional temperature 43, velocity perturbation <D, 

and an auxiliary function k characterizing the rotationality of the 
flow. With the use of this theory the action of heat sources on the 
flow is studied, and the heat delivery in a compressible flow at subsonic 
and supersonic speeds is calculated. The results show the effect of 
compressibility and the nonlinear cooling. Applications of the results 
to hot-wire anemometry are discussed. 


INTRODUCTION 


In a two-dimensional stationary field, consider a certain distrib- 
ution of heat sources placed in a horizontal stream of main velocity U, 
which can be subsonic or supersonic. The following points will be 
studied theoretically: 

(1) Action of heat addition on the motion of a compressible fluid 

(2) Effect of compressibility on the heat delivery 

(3) Nonlinearity in the relation between heat delivery and temper- 

ature difference 

For the investigations of these properties, the following problems will 
enter, and their analysis will form the main part of the present paper: 

(l) Development of a generalized- potential theory for the study of 
a compressible flow which is rotational and nonadiabatic (i.e., with 
heat addition) 
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(2) Temperature distribution about a heat source in a compressible 

flow 

(3) Distribution of the flux of heat delivered by the heat sources 

( 4) Perturbation of the flow by the heat sources 

This theoretical study will, on the one hand, illustrate the charac- 
teristic phenomena in a compressible flow with heat addition, and, on 
the other hand, it will find applications to hot-wire anemometry for the 
measurement of velocity in a compressible flow at subsonic and supersonic 
speeds. In order to simplify the calculations, without losing the 
essentiality of the problems, it will be assumed that the heat sources 
are distributed along a flat plate, which corresponds to an infinitely 
thin metallic ribbon, edgewise to the flow, so that no appreciable per- 
turbations due to the shape of the body will be present, and that pertur- 
bations are due only to the heat sources.-*- Actually a cylindrical wire 
customarily is used in hot-wire anemometry for practical reasons, but 
from the present calculations it will be seen that the form of a flat 
plate is preferable to the form of a blunt body like a cylinder in the 
hot-wire anemometry of a compressible flow, because the variations of 
the Mach number with the cooling are single-valued for the former body 
and double -valued for the latter body. After the calculations of the 
heat delivery by heat sources distributed along a flat plate, the heat 
delivery by heat sources distributed along a blunt body will be discussed 
also. In the present calculations the variations of the physical con- 
stants, such as the coefficients of specific heat and the coefficient of 
heat conduction, as well as the effect of viscosity and radiation, will 
be left out of consideration. 

This work was done at the National Bureau of Standards under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. 


-*-Since the viscosity is neglected, the terminology "plate," which 
occurs frequently in this paper, has the only meaning of indicating 
the location of the heat sources, without offering any perturbation to 
the flow, for example, in the form of a boundary layer. 
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GENERALIZED POTENTIAL THEORY FOR ROTATIOML MOTION OF A 
COMPRESSIBLE FLOW WITH VARIABLE TOTAL ENERGY 
Fundamental Equations 

By •writing the equations for the motion of the fluid, the following 
usual symbols will be introduced: 


x, y 


coordinates of a variable point in two-dimensional field 


ds = \/dx 2 + dy 2 
w vector velocity 

u, v components of vector velocity 

p pressure 

density 

absolute temperature 


R 

c 

X 

M 


ratio of specific heat for constant pressure c p and for 

constant volume c v (k = 1.4 is taken for air in numerical 
calculations) 

gas constant ^c p - c v = Cp — ~ ^ 

velocity of sound at temperature T ^\/kRT ) 
coefficient of heat conduction 
local Mach number (w/c) 

The following functions will also be introduced: 


P = 1 + m 2 


a 


= {l -\ 2 

1 


P = V* 


k-1 



k 
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2h - P 0 c p/^ 


H' = P^h 


2-k 

g = Ph'U = Poo k_1 hU 

The above symbols are frequently used in the present calculations. 
Other ones will be introduced wherever needed. The subscript o indi- 
cates the state reduced adiabatica l 1 y to stagnation, and the subscript oo 
indicates the state in the free stream. 

The total temperature at a variable point is given by the definition 




T + 



= T 0 + *(x,y) 


( 1 ) 


where T Q is a constant called isentropic stagnation temperature and 
■0(x,y) is the additional temperature due to an introduction of heat. 

For the sake of convenience the relations between the different temper- 
ature symbols are shown as follows: 


T t /T 0 = 1 + */T 0 = t 
T- fc/T = P 


( 2 ) 


Also, 


T/T 0 = Tp _1 


w c 


2c p T o 


= t(i - p-l) 


(3) 


In the present problem there are four dependent variables w, p, 
p, and i3. For these four variables four equations will be written by 
expressing the following four laws for the motion of the fluid: 

(l) The law of conservation of momentum 
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(2) The law of conservation of mass 

(3) The law of conservation of energy 

(4) The law of state of a perfect gas 

The conservation of momentum is expressed by the equation of 
momentum. In vectorial form it can be written as follows: 


l/2 grad w^ + rot w X w = grad p 


( 4 ) 


In a two-dimensional problem this equation for w can be separated into 


u 


u 


v as follows : 


du du 1 dp 

— + v — — 

dx dy p dx 

(*a) 

dv dv 1 dp 

dx dy p dy 

(4b) 


The conservation of mass of the moving fluid is expressed by the 
equation of continuity 


w dp 

div w = — — 
P ds 


(5) 


where ds = \/dx 2 + dy2 is the length of a segment of the stream l i n e. 

The conservation of energy is expressed by the equation of energy, 
sometimes called the equation of heat conduction 


pw 


T 


grad ^c v T + ^ w 2 j = - div (pw) + X div grad 
This equation can easily be transformed into the following one: 
pw • grad ^CpT + i w 2 j = X div grad T 


( 6 ) 
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In equation (6) T can be replaced by t from defi n ition (l), and 
the following equation can be expressed: 

1 

P = PoCT/To)^ 1 Si 


with 


Q (S 0 -S)/B 

S]_ = e 


where S is given by the definition of entropy 


T dS 


C P 



Then equation (6) becomes 


1 

div grad (t/p) = 2hSi(T/p) k_1 w * grad T (6a) 


where 2h = p Q c^/X. 

Finally the equation of state of a perfect gas is 

p = RpT (7) 

These five scalar equations (equations (4a), (4b), (5), (6a), and 
(7)) determine the five variables in the problem u, v, p, p, and ^ . 

The problem is now to reduce the number of equations by eliminating 
cert ain variables. Therefore a system of three equations will be dis- 
cussed which will determine the three principal variables u, v, and ■fl . 

The first equation can be obtained from the equation of continuity 
(5)> hy eliminating p and p with the use of equations (4a), (4b), 
and (7) as follows: 


div w = 


^ w£ 5w + w 5T 
c 2 5s T 5s 
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or again, by replacing T by S from equation (l), 


div w = 


w^ bw | w bti 
c 2 cis T c3s 


In the x, y coordinates this equation can be written in its 
expanded form as follows: 


u 2 \ du uv/ov Su\ L 

C 2 jbx ' c2\Bx by) [ 


v£\8v _ 1 / 8d \ 

c2 )by Ty bx by) 


( 8 ) 


Herewith the first equation of a system of three equations is obtained 
which will be discussed in the following paragraphs. The firBt equation 
will be called the equation of motion of the fluid. The second equation 
is the equation of energy already given by equation (6), and the third 
equation can be derived as follows from the condition of the rotationality 
of the flow which is given by the equation of momentum (4) . Simple cal- 
culations of the rotational term in equation ( 4) as a function of d 
and T-fc from equation (l) lead respectively to 


rot w X w 


1 grad log ( p k / p 




k - 1 


grad d 


(9) 


rot w x w 
Z2 


k - 1 


grad logp 


(p k /p) - 2kk ir loge Ti=1 ^ 2 ^ 


where grad log e 
expressed by 



is nothing else than the increase of entropy S, 


(p 1/k \ 

grad logg l— — J = grad (S/c p ) 


as can easily be proved from the definition of the entropy. 
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From equation (9a) it is seen that an isentropic (S = Constant) and 
adiabatic flow (T-fc = Constant) has 


rot w x w 

c2 


= 0 


Furthermore, by taking the rot of equation (9a) , making use of 
equations (2), and observing that rot grad A = 0, there is obtained 
finally 


. rot vxw , w^ q /2 

rot 5 = - grad X grad log e ^ 1 


( 10 ) 


This is the equation of rotationality of the flow. 

The two equations ( 8 ) and (10) determine the u, v field as a 
function of the temperature which, in its turn, is determined by the 
equation of energy ( 6 ). 


Generalized Potential Theory for a Rotational Flow 

A system of coordinates formed by the streamlines and their ortho- 
gonal trajectories is taken. The stream function is used for the 
streamlines, giving 


Po 

P 8y 

Po dt 

p 8x 


( 11 ) 


A generalized potential function cp and an auxiliary function k 
can be chosen such that 


U = K 

oy 


= K 

oy 


(12) 


V 
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or in vectorial form w = k gradcp. The function cp is called the 
" generalized" potential function, because it comprises the partic ular 
function for an " irrotational " potential flow by putting k = 1. The 
generalized potential motion can now be determined by three equations: 

The equation of motion in the form of a differential equation for cp, 
the equation of energy which gives the temperature distribution, and a 
third equation which determines the auxiliary function k. 

It is worth while to remark that the study of the distribution of 
a certain physical property (e.g., temperature distribution) about a 
body of more complicated form can be made possible by the introduction 
of such a set of functions cp and t in the generalized potential 
theory, because the use of such variables cp and ijr as new coordinates 
will transform the body into a flat plate. 

The third equation which determines the auxiliary f unc tion k can 
be constructed as follows. If the u, v field is considered momentarily 
as given, the function k will be determined by the following relations - 
derived by simple vectorial operations: 


rot w = -w x grad log e k 


(13) 


From equation (13) the rotational functions 
can be calculated as follows: 


rot w X w 
C 2 


and 


rot 


rot w X w 
C 2 


rot w X w 
C 2 



grad log e K + 


w • grad log e k 
c^ 


w 


(HO 


ro t r°t v x v . _ ^ 3>! x ^ w * . 


w • grad log e k 


w X grad log e K 


w X grad 


*w - grad log e k 


(15) 
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By comparing equation (15) with equation ( 10 ), the equation for k 
is obtained in the vectorial form as follows: 


w 2 w • grad log e k 

grad —x X grad log e k + 5 w x grad log- k + 

/w • grad log e k\ v2 - 1/0 . .. 

w X grad I — j = grad — X grad log e t 1 / (l6) 


The velocity potential 9 can be written in two parts aB follows: 

<P = ? + 0 (17) 

where cp is the perturbed velocity potential, q> = Ux is the velocity 
potential in the free stream, and 0 is the potential perturbation. 
Also, 


u = K (U + 8$>/8x) = U + u' 
u' = k 8<I>/8x + ( K - l)U >■ 
v = v' = k d<b/by 


(18) 


where U is the velocity in the free stream. 

After substitution of relations (l8), equation (8) becomes 


with 


L u 2 \a 2 * 

9 UV ^ , 

L v 2 ^ 2 * 

\ C 2/8x2 

CZ. » T 

c^ 8x8y 

l ‘ c2jSjr2 


f = 


M- 8t 

— w — 

KT 8s 





r (19) 


where 


ds = \jdxP- + dy 2 
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It may be useful to group the main equations as follows: 

(a) Equation of motion (see equation ( 19 ) ) : 


(± u2\ B% p uv B 2 $ / _ 

\ c 2 )bx 2 c 2 BxBy \ c 2 /By 2 


where f is given in equation (19) 

(b) Equation of energy: 


1 

A(h- lT ) = 2h(T/|a) k-1 w ^ 


a 2 b 2 

where A = g-p + ^-p. This equation is the analytical form of 
equation ( 6a) . 

(c) Equation of rotationality (see equation (l6)): 

, w 2 _ = . w • grad log e K 

grad — X grad log e k + — w X grad log e k + 


w X grad 


w • grad loge K 


= grad ^ X grad log e t 1 / 2 


If s and n represent the unit vectors along a streamline and 
along a normal, respectively, the following equations can be written: 


w • grad log e K = w 


5 loge K 


w x grad log e k = w 


6 log e K 
Bn 


, „p ^ n BM 2 d log e k ^M 2 B loge k 

grad VF X grad log e k = ^ ^ — — 

Bs Bn Bn Bs 
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Using these formulas and assuming 


3 log e K 

5s 


« 


d lQ g e 

5n 


K 


( 20 ) 


there is obtained from equation (l6) 

K m ]fr ( 21 ) 

In an adiabatic flow (d = 0) , the expression (21) reduces to the well- 
known value K = 1. Assumption (20) expresses the condition that the 
gradient of the additional temperature (or additional energy) must be 
much smaller along the streamlines than along their orthogonal trajec- 
tories. This condition exists in problems of the introduction of a heat 
boundary (problems of combustion), of heat sources, or of a hot body in 
a stream moving at a high speed. In all those problems the motion of the 
fluid is extraneous to the heat introduction. Otherwise condition (20) 
will not be valid, for example, in the case of an explosion, where the 
predominant motion of the fluid is created by the heat emission itself. 

In the present paper the problem of the rotational and nonadiabatic 
flow will be treated under the assumption (20) and hence with the value 
of k given by expression (21). Substituting expression (21) Into 
equation (19), the equation of motion becomes 

/ _ u 2 \d 2 <P _ uv; d 2 ^ /, v£\5^$ _ r 

\ c^dx 2 c 2 5x5y \ c 2 /5y 2 



Small Perturbations 

As a first approximation, equations (6a) and (22) can be transformed 
into linear ones by introducing the assumption of small perturbations; 
that is, it is aBBumed that the perturbations u'/u, v'/U, <&/(p, and 

■g/Tq and their differentials are smal 1 compared with unity, such that 
terms of second order are negligible. 
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Under -this assumption, the equation of energy (6a) and the equation 
of motion (22) become, respectively, 

(a) Equation of energy: 


where 


a 2 # a 2 # _ &9 

Sx 2 Sy 2 S Bx 


2-k 

"k-1 

g = Pco hU 


(23) 


(h) Equation of motion: 


2 a 2 $ a 2 $ 

ax 2 ay 2 = 


where 


(24 


1/2 


O' = (l - M» 2 ) 


f = i(l + kMn 2 )^- 

2 v ” 't d a x 


DISTRIBUTION OF TEMPERATURE AND HEAT FLUX FOR 
A HEATED BODY IN A COMPRESSIBLE FLOW 
Temperature Distribution 

The convection of heat is governed by equation (23)- It has the 
following solution: 


•3 = Ae gx K 0 (gr) 


(25) 



Ik 


MCA TN 2436 


where r = l/x^ + y2 and Kq is a Bessel function. For larger values 
of gr, Kq can he expanded into a series which will be given later in 

this section. The constant A can be determined if the heat delivery Q 
is known. Consider first the case of a line source. The formula of 
temperature distribution about a finite body can be constructed by 
considering a continuous distribution of heat sources on the body, and 
will be studied in the following section. 

The heat delivery Q should be equal to the total heat flux around 
a closed circuit formed, for example, by a circle of radius r enclosing 
the line source. The following formula can be written: 


cos 


where ds is a segment of arc on the circle and r and q are respec- 
tively the radius and the angle in the polar coordinates. 

The question is in what sector of the circuit the heat flow will be 
effective. At large velocities, the major part of the heat flow occutb 
through a very narrow sector so that the magnitude Q will be at first 
approximation independent of the limits of integration, provided these 
are taken of sufficiently large range. Therefore as a first approxi- 
mation there can be written in place of formula ( 2 6 ), for all kinds of 
velocities (subsonic and supersonic). 



Q 


= / as (- x s 


pCptFS 


Q = 2r\ J dTj^-^ + 2gp a f 1 -a cos rjj 


Putting 


low =1 


-if 


I-i(z) = 5 / dT ie 
J 0 


Z COS T) 


( 26 a) 


( 27 ) 


COS T) 



MCA TN 2^36 


formula ( 26 a) can "be written as follows: 


Q = 2AAgr 


-K].(gr) f 

J 0 


- gr cos Ti , 
dT^e 0 1 + 


- l)Ko(gr) J dne® 1, C0S 11 cos t] 


= 2itXAgr \ 


-K^gr^igr) + 


- 2(l - 


By the use of the following expansions. 


Kq(z) 


% e 


Vl'l 1 - 8?) 


KlW - £ « --Vlft ♦ £) 


Io(z) " S 1 + s) 


Il(z) = f» -7==(l ~ 

1 dz ^/2^\ 8z/ 


formula ( 28 ) becomes 

Q = 2jt\APoo 

Herefrom A is obtained, and from formula (25) there is obtained 


-1 


4 = A eSXK o(^° 
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Flux Distribution and Heat Delivery 

The distribution of temperature about a line source is given by 
formula (29). This formula is valid for subsonic and supersonic flows. 
Consider now a flat plate of length 4a. Let q(£)d£ represent the heat 
flux from an elementary segment d| of the plate. The temperature 
contribution by this segment at a variable point x,y is 

dS = — dgq(6)e S ^ X "^K 0 ^g|/z 2 + y 2 ) 


where z = x - 5- By integrating. 


d(x,y) 


M- CO 

2jc\ 


J dgq(l)e S ^ X- ^K 0 (g\/z 2 + y 2 ) 


(30) 


Formula (30) gives the distribution of temperature in the flow region 
about a plate. This formula contains a function q which will be deter- 
mined, in the calculations which follow formula ( 30 ), by the integral 
equation ( 31 ) for the surface condition on the plate. 

On the surface of the plate there can be written for 0 £ x S 4a 


^ = / d|q(|)e g(X ^KqJ^x - lT] + 

0 


2flX 


f d|q(i)e" s(l_x) K 0 [^(i 


x) 


(31) 


or, by using the asymptotic expansion of 


K^z) for large values of z. 


w 




e -2g( |-x) 


n 


(31a) 



MCA. TN 2436 


IT 


The second integral on the right-hand side decreases rapidly -with £, 
and therefore as a first approximation it can "be written as follows: 



d£q( |)e 


-2g( |-x) 


it 


|2g(5 - xi 






For the sake of simplification of writing, introduce 


where 


qp(x) = 


1 . 
g* 2*1 q(x) 


£* = ^2g/rt 




The relation (31a) between the heat flux and the wall temperature then 
becomes 


or 


"®w - 






qpU) 





qid) 

)[* - 1 


(32) 


First the integral equation (32) will be solved for the general case of 
a variable wall temperature Later the particular case will be 

derived of a constant \ which is true, if the plate is supposed to 
be a very good conductor- 
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The integral equation will he reduced to a linear differential 
equation of first order by a differentiation with respect to x. First, 

put x - | = £ 2 in order to prevent the singularity which would occur 
■when | = x. The differential with respect to x of the integral in 
the right-hand member of equation ( 32 ) is 



q 1 (x - £ 2 ) 


+ 2 



d£qi(x - £ 2 ) 



4p( i) 
Vx - i 


where q-j^C l) = — 

By applying this rule to the differentiation of equation (32), there 
is obtained 


q-^x) 




( 33 ) 


This formula is fulfilled by the flux at the distance x. The flux 
at the distance Xg will fulfill 


q x (X2) + g* 


qp(o) 

^2 


+ e. 




L 


*2 


dl 


^(6) 
)p2 ~ 1 


= Vw( X 2) 


( 3*0 
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"“l/2 

Multiply every term of equation (3*0 "by dxgg^x - * 2 ) ' > 

integrate with respect to X£ between the limits 0 and x, and sub- 
tract from equation (33)j it follows: 



(35) 


The first integral on the left-hand side of this equation has the 
value it. In the double integral the order of the integration can be 
inverted as follows: 



The last integral with respect to Xg has the value n. Hence the 
value of the double integral becomes 
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Substituting this into equation ( 35 ) 3 


q-j_(x) - 2gq 1 (x) = F(x) 


where 


F(x) = - g# q 1 ( 0 )x- 1 / 2 + g ^ w - g 2 f ^ 

'Jr, 


K(*2) 


The integral equation ( 32 ) is thus transformed into the linear 
differential equation ( 36 ). Its solution is 

q-jU) = qi (0)e 2gX + J d|F(i)e 2s(x_l) 

= q 1 (0)e 2gx ^l - erf \/2gx) + q^x) ( 


where 


2 pi 2S X _x2 

erf \j2gx = pj=( J tee 


and q^' (x) is the flux due to the gradient of the wall temperature 


qi'(x) = I" d£F'(|)e 
^0 


2g(x-i) 


with 


F'd) - s*«„ - g. 


r ^ 

J o 


^w(*2) 
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In the particular case of a plate which is a good conductor, the 
wall temperature can he assumed constant; then formula (37) becomes 
simply 


ll(x) = q-^Oje 2 ^! - erf )j2gx'j (37a) 

or, by returning to q, 

q(x) = q(o)e 2gx ^l - erf ^2gxj 


This relation is plotted in figure 1, where \j2gx is denoted by u, 
which can be expressed as 


2-k 

"Stay nr 

u - v 5 ^ te 


where P is the Peclet number as will be defined later. 

The constant q 1 (0) in formula (37a) can be determined by putting 
x = 0 in equation ( 32 ) which gives: 


*L (0) ‘ 2irXg* « (0) * 


(37b) 


Now pass over to the calculation of the heat delivery. By definition. 


or 


% = P dxq-^x) = f dxq 1 (0)e 2gX fl - erf ]/2gxj 

Jo Jo 

Ql/ q l(°) = dxe 2SX ~ /** dxe 2gX erf \j2gx 

Jq Jq 
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The first integral of the right-hand member is evident. By changing 

the variables = 2gx, the second integral of the right-hand member 
can be written as follows: 



Hence , 


Q % 

qCoT q^ToJ 


2g 







where, for ga » 1, ct^ is 

a-, = ~ ]/8ga 


(39a) 


Define a Nusselt number N such that 


Q = NAA 


where A is the area per unit length 8a. Hence, 

Q 


N = 


(4o) 
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From expression (37t>)> 


q( 0) = 2\<S V jjfi 

Hence, 

N = 2^; = ff iA» = f 


(41) 


(42) 


The above calculations give a picture of the distribution of the 
heat flux. However, in the order of approximation adopted, it seems 
that the pattern of distribution of the heat flux does not influence 
very much the heat delivery. This is readily seen, since the expression 
for the heat delivery (formula (42)) can alBO be obtained by the 
following more approximate method which does not take into account the 
details of the flux distribution. To this end, if in equation (31) 
a mean heat flux q is introduced and x is put equal to 4a, 
the following equation can be written: 


' s ^SS?/ lta «= g “ ta ' l)K olj^- »] 




>4ga 


2 nX 


qg 


J due u K 0 (u) 


(42a) 


This integral is known and has the value \j8nga - 1. Then, 


q = 2X,T3 w JTgp co " 1 (\j8nga - 1^ _1 
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and N which is defined by 


4aq 

2 ?^ 


becomes 



By comparing formulas ( k2) and ( 43) , it is seen that formula (43) gives 
a very good approximation to formula (42) which was obtained from the 
considerations of the exact distribution of heat flux. 

A generalized potential theory has been developed and three differ- 
ential equations have been obtained which determine the three variables: 
The distribution of temperature difference i3, the potential pertur- 
bation <&, and the auxiliary function of rotationality K. As a first 
approximation, the equation has been linearized by neglecting terms of 
the second order in perturbations u'/U, i>/cp, and i^/Tq (see section 

Small Perturbations). Then the equations can be solved for t 3 (section 
Temperature Distribution), the heat flux can be calculated (this section), 
and hence the heat delivery can be obtained as in formula (4-3). How- 
ever, the results show only a linear relation between the heat delivery 
and the temperature difference. In order to show the nonlinear char- 
acter, a second approximation must be pursued by keeping the second- 
order terms of perturbations in the calculations of heat delivery. To 
this end, the distribution of temperature difference $ obtained in 
the first approximation will be used to calculate first the potential 
perturbation <I> (two following sections), and herefrom the velocity per- 
turbation u'/U- An "effective velocity" which is the average perturbed 
velocity of the fluid along the plate can be introduced in the expression 
of the heat delivery which will then be of nonlinear nature (see first 
section under Nonlinear Cooling of a Heated Body at Supersonic Speed) . 

It is to be remarked that the process mentioned above will give a 
second approximation with sufficient accuracy without entering into the 
integration of the nonlinear partial differential equation in the form 
of equation (6a). 
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ACTION OF HEAT ADDITION ON MOTION 
OF A COMPRESSIBLE FLOW 
Flow Perturbations by Heat Sources 
at Subsonic Speed 


In this section the flow perturbation at a subsonic speed by heat 
sources distributed along a flat plate will be investigated by using 
equation ( 2k) . By means of the transformation 



m 


equation {2k) can be brought to 


a Poisson' s equation of the form 


d 2 ^ d 2 ^ 

dxj 2 dy ^ 2 


fVa " 2 


(* 5 ) 


For the sake of mathematical convenience this equation for the 
motion of the fluid about a plate can again be transformed by conformal 
representation into an equation for the motion of the fluid about a 
cylinder of radius a, using the formula of conformal transformation 



where z = + iy^ represents the plane of the flow about the plate, 

and Z = X + iY the plane of the flow about the cylinder. By trans- 
forming equation (lj-5) into an equation for a cylinder, there is obtained 
the following equation: 


d 2 ^ B 2 ^ 

"ix 2 " + ai 2 " 


fVa" 2 




{k6) 


v 
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TM 3 equation can be made dimensionless by dividing X and Y by 
a, an d <5^ by 2Ua. When R and T} are taken as the polar coordinates 

' of the Z -plane, and 


dz 

dZ 


= 1 - 2R“ 2 cos 2x\ + R“ 


the right-band member of equation (46) becomes 

f ± = i va _2 (l + kM^l - 2R -2 cos 2ii + 


(47) 


where dd/dx should be obtained by differentiating expression (30) for 
the temperature distribution. However, the latter expression's being 
of integral form w ild make the Poisson's equation (46) rather cumbersome 
to handl e. As a convenient but not as a necessary step, the perturbation 
function <tj_ is considered as a superposition of elementary perturba- 
tions corresponding to elementary heat sources. First calculate the 
elementary perturbation produced by an elementary heat source from 
equation (46). An integration according to the proper distribution of 
heat sources will give us the resultant perturbation. 


Following the above procedure, consider first a heat source of 
intensity Q, placed at the origin of the coordinate (center of the 
plate); then dd/dx can be calculated from formula (29), and expreB 
sion (47) can be written as follows: 


f ± «* B-jV 


)Jsax± 


-gar i (l-xi/ri) 

; (1 - xp/ri) 


dz 

dZ 


(47a) 


where 


B 1 = 


4^2jt 


(l + kM 2 V -2 

\ 00 / 00 T 


a" 2 (l + 


Q 

2Xd. 


ga 


w 


It is worth while to transform this formula into the polar coor- 
dinates R and r\ of the Z-plane. To this end, note first that the 
exponential term in the function f^ can be reduced to the form 

-constant X ga(R-l) 2 

fp ~ e 



MCA TN 24-36 


27 


Since large values of ga must be dealt with, the function f^ keeps 
a significant value only for £ = R - 1 « 1. Therefore the approxi- 
mation £ « 1 can be introduced in the transformation process. Under 
such an approximation, the variables x^ and y^ in terms of £ and 
tj become 


x^ = 2 cos T) 


with 


y i = 2£ sin t] 

r i = {h 2 + y ± 2 


It follows that, at a variable point downstream of the source for 

-rt/2 « t] « + it/2, 


V2g< 


;a cos 


(V7b) 


Now the integration of Poisson's equation (46) for the Z-plane will 
be considered. It is to be remarked that the effect of the presence of 
a right-hand member in equation (46) is equivalent to that of a certain 
continuous distribution of sources throughout the region of flow, the 
intensity dl of the source distribution in an element of area l d l dX 
being 


dl = l dZ dX 
2 « 


where l and x are the polar coordinates. The potential function 
induced by this distribution of sources can be calculated by the method 
of image (references 1 to 3) • 

Consider first a unit source placed at a point Zq = Ze ix in the 
flow field outside the cylinder. The velocity potential induced by this 

source at a variable point Z = Re^ is equivalent to the velocity 
potential of the unit source placed at Zq, its unit image placed at the 

*1 4 V 

point Z r = reciprocal to Zq, and a unit source of intensity -1 

placed at the center of the cylinder. When the potential function and 
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stream function of the perturbation induced by the uni t source are 
denoted by and $*, respectively, there can be -written 


® + ilf* — log e ^ Z - Zq^ + logg ^Z - Zj>^ - logs Z 

By ta k in g the real part of the above complex equation, there is obtained 

= | log e (z 2 - 2ZR cos (X - Tj) + R 2 ) + 


||l - 2(ZR) -1 cos ( X - tj ) + (ZR)" 2 ] 


Evidently this expression can be written in the two following identical 
forms: 


= § log e 7 2 + \ loge [l ~ cos (x - T]) + + 

~ log e Jl - 2| 2 cos (X - tj) + 

= § loge R 2 + \ log e Jl ~ 2 ^l' cos ( x - B) + (Si 1 ) 2 ] + 

§ loge [l - 2 |g cos (X - Tj) + £g^| ( 48 ) 


where, for the sake of simplification, 

ll = R/Z 

h' = Z /R 
^2 = (ZR )" 1 
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The total induced velocity potential can he obtained by integrating 
$* for the whole region outside the cylinder as follows: 



00 2jt 

J d2 J dXZf ± ( l > X)$*( Z, X) 


m 


Thus the solution of Poisson equation (^6) is reduced, to the solution 
of the double integral {^9). 

Substituting expression ( - 48 ) into formula (b9) , the double integral 
can be written as follows: 


®i = | -J dl J dXZfijl log e R 2 + | logg [l - 26i * cos (X - ij) + 

+ | log e [l 2^ cos (X - tj) + Ig 2 ] 

rv 00 n2^ 

| -J dl J dXZfJ| log e * 2 + | log e f 1 " cos (X - tj) 

R 0 


ll 2 ] + | log e [l - 2 1 2 cos (X - t]) + 5 2 2 J 


(50) 


Only the velocity perturbation at the surface of the cylinder (R = l) 
is of interest; this is 


l + 


[l - 2| cos (X - tj) + g 2 ] (50a) 


with |=1 /l. 

It can be proved that for a flat plate the perturbation upstream 
of the heat source decreases very rapidly with the distance, in the 


*io “ h f dl f 51 ^1 lo Se 
J 1 J 0 

/ °° n2rt 

anj axf ± loge 
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proportion e I and hence it is negligible with respect to the 
perturbation downstream. This means that the variations of f^(Z,X) 

have to be restricted to the half plane corresponding to -it/2 and 
ic/2. Therefore the limits of integration in formula (50a) will be 
-it/ 2 and it/2. In order to carry out the integrations of the right- 
hand member of formula (50a), the following simplifications can be 
introduced: 

(1) The variable Z can be changed into 1 + £ with £ « 1 

(2) Since log e Z ** £, the first double integral of the right-hand 

member of formula (50a) can be neglected in comparison with the second 
one 


(3) The order of integrations in the second double integral of the 
right-hand member can be interchanged; by first integrating f^ with 
respect to 


d£f ± (Z,X) = B-.V (ga)~ 2 tgX 

± 2 \J2 


Hence formula (50a) is reduced to 


® io * B]V (ga )" 2 dX log {2I1 - cos (X - tj )1 

-it/2 l L J 


tgX 


= BjV r (ga ) -2 r 
J 0 


n /2 sin 2 X 1 1 - cos (X - ti) 

dX logg 


cos X sin X 1 - cos (X + tj) 


As the logarithmic term has predominant absolute values in the 
proximity of X ~ tj, there can be written by approximation 


“ BlV ip (ga) 


0 

-2 sin n 


COS T] Oq 


1 


,n/2 


dX 


sin X 


log. 


1 - cos (X - T|) 
1 - COS (X + T|) 
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Now the integral on 

the 

(1) For 

H = 0 

it 

(2) For 

'■I 

it 


right-hand side has the folio-wing limit values: 
is zero 
is 



dX 


1 

sin X 


log e 


1 - sin X 
1 + sin X 



(See reference 4.) For intermediary values of T], the approximate 
n 2 

form - — sin T] is suggested. Hence, 


io 


Rj -B-|V 


/ s-2 sin^ri 

— (ga) L 

8v2 cos t) 


(5D 


The distribution of the potential perturbation $> 0 on the compres- 
sible field along a plate of length 4a can be obtained by multiplying 
formula (5l) by 2Ua, by using relations (4-4), and by replacing 
cos T] by x/2a; then. 


4> 0 = -Bp ^j| (ga)~ 2 Ua ^|l - (x^a) 2 ] 3 ^ 2 
or, by substituting for Bj, 

*0 - -f B 2 IKS- 1 JL- f |1 - (x/2a) 2 ] 3/2 (52) 

where Bg represents for the sake of simplification 

B 2 = ^ a" 2 (l + kM ro 2 ) ^ (53) 


Formula (52) gives the potential perturbation by a single heat source at 
a point x from the source. 
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Now consider a continuous distribution of heat sources such that 
the intensity in the region between | and I + d| from the leading 
edge of the plate is q( I) d£: Thus the elementary perturbations by 

q( |) d£ along the downstream part of the segment of plate which has a 
length 4a - £ can be obta i ned from formula ( 52 ) by replacing 2a by 
4a - | and x by x' - £, where x' is the distance from the leading 
edge. For x' - | ^ 0, 


<a<K 0 ( t) 


it _ 

= B 2^oog 


q(S) dg 4a - g 

2Xd w x 1 - | 


1 



The total downstream perturbation corresponding to a continuous 
distribution of heat sources along the whole plate is obtained by an 
integration of the above expression with respect to | between the 
limits 0 and x’ as follows: 


9 


o 


y B 2 U — — 

4 2 2g 24^ 





(54) 


In order to calculate the integral in the right-hand member of the 
above expression, the following simplifications are introduced: 

(l) The asymptotic form of the formula (37a) for flux distribution 
will be used 


q(l) w q(o) 


^2 Jtg | 


where 


/ \ A 1 CO 1 V 1 
q(0) ~ ( 2g 24^ 


(2) A change of variables ^ 1 = cos u transforms expression (54) 

into the following type: 

a, - sln * - r(„) 

a “It 


COS •‘r 




cos u 
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where F(t|) is a trigonometrical function which appears in the course of 
the transformation and v is a parameter which turns out to he 
v = x'/^ a ' The detailed form of F(ti) will not he written out here. 

As the fractional term has predominant values in the proximity of 
‘ cos u v, there can he written hy approximation 


iit/2 

"cos~^v 


/ 3,1 


sin 


\Jv - cos u 


F(ti) k F(v 


jt/2 


>/: 

cos 


■^V 


dTi 


sin tj 


Vv - cos u 


The integral on the right-hand side can he immediately calculated, and 
there is obtained 


$o = -j® 2 y a ~ v 1/ ^ 2 ( 1 ~ v )(l + v) 3 / 2 

The gradient at the center of the plate will he taken as the effective 
gradient e which is 


c _ 1 d$ o 
U dx 


3€ 


\fiBr 


a 


= ^-\/ita _2 (l + kM^ 2 ) — —= 

T Q \jQga 


( 55 ) 


It is seen that in a subsonic field the fluid increases its velocity hy 
passing over the heat sources. 


Flow Perturbations hy Heat Sources at Supersonic Speed 

In this section the flow perturbations at a supersonic speed hy 
heat sources distributed along a flat plate will he investigated hy using 
again equation (2^). The equation will he first solved for a single 
heat source placed at the origin. The perturbations of the flow about 
a flat plate will he derived hy a superposition of heat sources as done 
in the preceding section. Equation (2b) can he made dimensionless hy 
dividing the coordinates hy a, and the potential perturbation hy 2Ua, 
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where 4a is the length of the plate and U is the free-stream velocity. 
Thus, 


where 


_ 1 _ _ d 2 g> 

cr 2 dx 2 dy 2 


= -f 




a = (M ^ 2 - l ) -1 / 2 = 


■fcgric 


t ] 0 = Mach angle 


(56) 


By means of a change of variables 


= crx 

y± = y 


(57) 


equation ( 56 ) can be transformed into the following form: 


5 2 $ _ f 

9y ± 2 SxjL 2 


( 56 a) 


Equation ( 56 a) has its characteristic lines given by yj t = C; these 

are two families of straight lines making an angle of 4-5° with the 
x^-axis. 


The problem is to find a solution of equation ( 56 a) such that it 
satisfies the following two boundary conditions: 

(a) The term b$/<5y± vanishes at y^ = 0 

(b) The perturbation $ vanishes along the front shock waves, which 
are formed by the characteristic lines originating at the leading edge. 
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In particular, tlie solution at various distances x^ on the x^axis is 
of special interest. 

The method of characteristics will he used for the integration of 
equation (56a). Therefore, take a point x-^ on the x^-axis, and draw 
a line x^p parallel to Otj * (see fig. 2 ). In the plane x^y^ apply 
the Gaussian theorem in the form 

JJ ^ - ll) - f (D + v ^ ( 58) 

Oxpp ' 

where the double integration is extended over the area Ox^p, and the 

line integral along the contour of this area in the positive direction. 
If in equation (58) there are inserted U = y^ and V = c^/cbc^, it 
follows from equations (56a) and (58) that 

+ ‘-fl ^ ^ (59) 

Ox^p 

It is evident that dx.j_ = dy^ on the line Op, and dx^ = -dy^ on the 
line X-jp. Consequently, 
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and it follows from equation (59) that 


^XjL - 24> p + $ 0 + J 1 dx ± ~ = - JJ dx ± dy-jf 

0 Ox_j_p 


Nowy from the boundary conditions. 


' da \ 

fril 


= 0 


yi=o 


®p = ®o = 0 


Hence, 


$X± = ~ II ^ dyif 

Ox ± p 


(60) 


The above formula can be transformed from the x^yp -plane into the 
x,y-plane according to relations (57)* Omitting the details, the result 
is 


®x 



ayif(x iy i) 


(61) 


where l/a is the ratio between the two areas of integration in the 
respective planes 


1 _ 6(x,y) 
a S(xi,yi) 



MCA TN 2^36 


37 


For a heat source placed at the origin, f is given by expres- 
sions ( 56 ) and ( 29 ), and there results 


f ft* B- 


, 1 -gar(l-x/r) 

\I&BX 


(1 - x/r) 


with 


V 


h\[2rt 


(1 + kM* 2 )^ 2Xd 7T ga 


o w 


( 62 ) 


It is to be noted that an integration of f with respect to y^ 
means an integration with respect to 


v 


^gaCr-L - x) 


while keeping x constant. After some simplifications, 

T X 1 , B-,' , pV ga(\/2-l)xi 

/ ^l^Xl^l) = — ~ (ga) — / dvr 2 e'^ 

J 0 ft X 1 4 


By substituting this into formula (6l), it follows that 


0 r*/ 2 dx, n\/ga(\/2-l)x 1 « 

*x=-!vw- ! J ^/ v 


and, introducing a new variable u such that 


u = \Jga.[\[2 - l)x x 
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there results 



( 63 ) 


where 


\J2 - 1 

*2 = S a § X 

By interchanging the orders of integration, there can he written 



2V2 


Bi'(ga)' 




2 -v^ 

dw e 


( 64 ) 


For very small values of <& x is small, hut it increases 

rapidly until a value of X 2 in the order of l/(ga) is reached. Then 
the integrals take practically their asymptotic values which are 


and 


J o k 



2 -v 2 
dw e 


log e v = Constant 


00 
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Hence formula (64) 


becomes 


. \Z* B l' . , 

<t x (ga) 

0 \]2 


-2 


l°g e yfe - Constant 


As the constant does not play a material part in the velocity pertur- 
bation, it can be disregarded, and there is obtained 




\Jn Bj' 
a \f2 


-2 


(ga) ~ log e ^jx , 2 


and, by substituting for Bj' from expression (62), 




x 


JL_ 

8g 


(1 + kM^JpJga)" 1 


2X'9, T l0ge \/*2 


T o 


(65) 


Now consider a continuous distribution of heat sources over the 
length 4a of the plate, such that the intensity of the source on the 
segment between | and § + d| is q(|) di, where q has the distrib- 
ution law given by formula (37a). The elementary potential pertur- 
bation d<I> x due to q(0 di at a point situated at a distance x from 

the leading edge on the x-axis can be easily obtained from formula (65) 
by changing in formula (65) Q into q(| ) d£ and x into x - £ . 

Then, 


d<t> x = - , 
x l6cr 


7- (1 + kMj^Jga) -1 tr log, 


2Xd. 


V2 - 1 / ^ 

ga o — ( x “ I) 


w 


If the perturbation upstream of the source is again negligible as 
in the subsonic case, the perturbation $ corresponding to a contin- 
uous distribution of the heat sources along the total length of the 
plate can be obtained by integrating the above expression between the 
limits 0 and x of the variable £. Thus, 
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to 


*x = 


1 60 


(1 + kJ^ o 2 )(i oo (ga) _1 — f* di 

T o '■'n 


fa E ^> loge 

\J2 - 1 , m 

ga (x - |) 


2 

— _ 


( 66 ) 


For the sake of simplification in the following calculations, the 
integral in formula (66) will he represented by 


f d| log e 

ga ~ — (x - i) 


2 


(66a) 


Generally in a supersonic flow ga is very large, so that q( |) 
decreases very rapidly with i; therefore expression (66a) can he written 
in the following approximate form: 


_ . / va - 1 \ r x 4(6) 

% = i°se (*» ~i- J o 

The integral of- the right-hand member has been calculated in the case with 
formula (to) and has the value 


^(5) 2 1 — 

di — — \/2gx 

2Xd w \/n ^co 

where x is now the distance in absolute value from the leading edge. 
Hence formula (66) becomes 



x 


-~p - (1 + ^Moo 2 ) — (ga) -1 \/2gx log (g 

8^a T q \ 


V2 


2 


x 


(67) 
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The perturbation $ x increases in absolute value from the leading edge 
to the trailing edge with an average gradient 


^x 

dx 


k\[n 


r (1 + M„ 2 ) — 


T„ 



m - 1 


4ga U 


(68) 


It is seen that in a supersonic flow the heat sources will produce a 
decrease in the velocity of the flow which passes over them. 


NONLINEAR COOLING OF A HEATED BODY AT SUPERSONIC SPEED 
Nonlinear Cooling of a Flat Plate 


The role of the heat sources along a body in a compressible flow can 
be I m agined to heat a thin layer of air which perturbs the motion of the 
main flow, just as if the body were deformed instead of being heated. 

The average velocity which passes along such a deformed body is then an 
"effective velocity" defined by 


or 


U 


eff 


= U + 


^x 

dx 


U S f f /U - 1 + § 




U dx 


where d^/dx is given by formula (55) in a subsonic flow, and by 

formula (68) in a supersonic flow. For the sake of convenience of 
writing, introduce 


Peff ~ U .l^c 
U U dx 


( 69 ) 


as the correction coefficient of the effective velocity 


Ueff/U = 1 + e 
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Introducing this effective velocity in the expression of heat 
delivery ( 43) , 


N = 



^8irga( 1 


+ e/2) + 1 


(70) 


where 


2-k 


k-1 

g = ^co 


hU 


In problems of heat transfer, the following parameters play an 
important part: 

(l) The Russelt number N as introduced previously in formula (4o) 
is defined by 


Q = NX A \/d 


where Q is the amount of heat delivered in unit time by the unit length 
of an immersed body across an area A, d w is the difference between the 
temperature of the body and the isentropic stagnation temperature, and 
d is a representative length. For a flat plate of length 4a, N is 


N = 


Q 

21 ^ 


(2) The Reynolds number R = p 0 Ud/r), where U is the free-stream 
velocity, p Q is the isentropic stagnation density, and i} is the 
coefficient of viscosity 

(3) The Rrandtl number Rr = 2hq/p 0 , where 2h = p Q Cp/X 

(4) The Peclet number P = R x Pr = 8hUa 

It is to be noted that if the Prandtl number is a constant, the 
Peclet number is proportional to the Reynolds .number. If the results 
are to be expressed in terms of P, the quantity 8ga which occurs 
frequently in this paper becomes 


8ga = Pp ro 


2-k 

k-1 
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Hence an alternative form of formula (70) for the heat delivery is 
obtained as follows: 

k 

2N = fp(l + | 2(k_1) + [i^ -1 (70a) 

where e is given by formulas (55) and (68) for the subsonic and the 
supersonic speeds, respectively, as follows: 

(a) In a subsonic speed: 


6 = 


o 2-fc 

MI nr 1 + KM « ,1 2(k-i) 

^ V 71 o ^00 

16 1 - Moo 1 2 


_1 

p 2 « 


Wj 


•/T 0 


(71) 


(b) In a supersonic speed: 


6 = (1 + kM 00 2)(M C0 2 - 1) n, 


1 2-k 

2 2(k-l) 1 

— log e 


\/rtP 


2-k 


>fa - 1 g(k-i) ^ 
— T— >*« V? 


(72) 


By substituting formulas ( 71) and (72), formula (70a) can be written 
also in the following form: 


2N 


= C Q ^P + Cp - (c 2 log e ft - C3 ) ^ 


(70b) 


where N is the Nusselt number, P is the Peclet number (it is equal 
to the Reynolds number, when the Prandtl number is unity), is the 

difference between the temperature of the body and the isentropic stag- 
nation temperature T 0 , and finally Cq, Cp, C 2 , and C3 are the 

following functions of the free-stream Mach number M (see figs. 3 
and 4) . 
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(a) In a subsonic flow: 


c 0 


k 

- 2(k-l) 


C]_ — l-ioo 

c 2 = 0 


> 



1 + kMpcp 

1 - K? 



(73) 


(b) In a supersonic flow: 



k 

- 2(k-l) 


C 1 = Hco 

C 2 = | (1 + kM^KM* 2 - l) 1 / 2 ^" 1 


C3 = C, 


2(k - 1) 


l°Se ^00 " 1o Se 


I 


\J2 - 1 


(73a) 


where p ro = 1 + k M^ 2 , and k is the ratio of the coefficients of 
specific heat. 

From the formula (70b) for the heat delivery of a plate, the 
following properties can be seen (figs. 5 and 6). 


(a) Subsonic flow: 

(l) IT increases with tS w . The rate of change does not depend 
on F. Since Q = 2M&W, Q does not vary linearly with \ (fig- 5) 
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(2) N decreases as increases 

(3) N varies linearly with \f¥ 

( k) For very moderate M^, even for ly^, w 0, formula (70b) 

reduces to 


2N = + 1 + ^ d v /T 0 

It reduces to the formula of King (reference 5 ) only if $ W /T 0 vanishes, 
(b) Supersonic flow: 

(1) N decreases as 8^ increases. The rate of change depends 

on log e \[F and therefore is larger in absolute value than the rate of 

change in a subsonic flow; Q does not vary linearly with t9 w (fig. 5) 

(2) H decreases as Moo increases. This is true for a flat 

plate with a sharp nose. It will be seen in the following section that 
for a body with a blunt nose different aspects will appear 

(3) N does not vary linearly with \fP, because the temperature 
term in formula (70b) contains a function log e \/¥. However, the devi- 
ation from the linearity is small owing to the small magnitude of \/^ 0 • 


Some Remarks on Heat Delivery by a Body With 
a Blunt Hose in Supersonic Flow 

In the preceding sections, a certain continuous distribution of 
heat sources immersed in a flow of velocity U has been considered. By 
the use of a generalized potential theory developed for a nonadiabatic 
and rotational fluid, the action of the heat sources on the flow has been 
studied, and the heat delivery has been calculated. 

In order to show in the theory the fundamental action of heat on 
the flow, and to obtain the characteristic physical features, it is 
believed that the extraneous action of the shape of the medium, that is, 
the metallic body introducing heat to the fluid, must be eliminated in 
a first study. Therefore, in the present calculations, a fine plate has 
been chosen which does not offer any shape perturbation to the flow. 

This circumstance approaches the theoretical distribution of heat sources 
without interference from an extraneous body. Fundamental investigations 
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of the physical phenomena under such a simplified circumstance may serve 
as a basis for the study of the phenomena under other more complicated 
circumstances, for example, the action of heat combined with the action 
of an extr ane ous body- Hot-wire anemometry in an incompressible fluid 
takes us ually a cylindrical wire for the introduction of heat to the 
fluid. No difficulty has ever arisen since its body perturbation is 
known by elementary theory. In a supersonic flow, the body perturbation 
is more difficult to calculate. However, for the study of the heat 
delivery an approximate estimation can be made by starting from the 
following considerations. 

A blunt body like a cylinder will produce a detached shock wave 
which is strong on the front region and weak on the two branches. Behind 
the shock wave, the flow about the cylinder is subsonic on the front 
region nnd supersonic on the rest of the region. From the present cal- 
culations (see section Flux Distribution and Heat Delivery) it is known 
that the air in the front region gives the largest contribution in the 
heat delivery. That part of the air, heated under a subsonic state, 
passes the sonic zone and leaves the cylinder with a supersonic velocity. 
While the detailed picture of the phenomena has to be studied more rig- 
orously, a rough estimation of the heat delivery can already be obtained 
by giving to that contributing part of the air an effective velocity 
after a Mach wave, and an effective subsonic Mach number Mp which can 
be derived from the unperturbed Mach number M^ before that part of 
the shock wave which is approximately normal. The known relation between 
such Mach numbers is (see fig. 4) 



1 + 


— Moo 2 


kM^ 2 - (k - l)/2 


(74) 


The temperature distribution about a cylinder is given by 


6 2 t 3 + 5^-8 2g 

dq> 2 ^2 U &p 


(75) 


This is the general form of equation (23). By choosing a system of 
coordinates 9 and ty, this distribution reduces to the temperature 
distribution in a unif orm stream flowing along a flat plate, as given 
previously by equation (23). The same method used in the sections 
Temperature Distribution and Flux Distribution and Heat Delivery can be 
applied here in order to calculate the temperature distribution and 
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finally the heat delivery for a cylinder. Omitting, the details, the 
result is 

Q = + 1 ) ( 76 ) 

The limits of variation of the potential function are 0 and 0 Q . In 
the case of an elliptic cylinder of semiaxes (a,b), 

f3 0 = 2(a + "b)U 


As P Q = If-Ua for both cases of a circular cylinder of radius a and of 

a plate of length ka, it is seen that the expressions of the heat 
delivery for both cases without heat perturbation must be identical. It 
is also to be noted that, for M-j_ = 0, the expression of the heat deliv- 
ery reduces to that found by King. Further, for a circular cylinder, 

N = — Q — 

c ttA,i3 w 
P c = 4hUa 


( p c)eff 


4hU 


eff 


a 


Hence formula ( 76 ) becomes 


N c - c 0\j n ( F c) e ff + (^l) 


-1 


( 77 ) 


The mass of air on the front part of the cylinder, heated at a 
subsonic state characterized by the subsonic Mach number M^ will leave 

the cylinder with a supersonic velocity which is approximately the 
effective velocity after a Mach wave 


u e ff = U(1 + e) 
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where e is given "by expression (72). By substituting expression (72) 
into formula ( 77) ^ 


Nc = fcopT + | Cl - ^(<=2 lo «e fc ' 


(78) 


where the C's are the following functions of the Mach numbers: 


Cq(Mi) 


k 

”2(k-i) 


Cl(Mi) = ^i" 1 


C 2 (M X ) 


c 2 (iO 


Cp(Ml) 

CqOO 


C 3 (M 1 ) = C 2 (Mi) 


2(k - 1) 


loge l^oo " l°g e 


\!2 


Formula ( 78 ) remains valid for a subsonic flow, provided the values 
of the C-f unctions are taken from expressions (73). 


The amount of heat Q does not vary with (see fig. 7) j although 

the deviation from linearity is not so important as for the case of a 
flat plate. 

For small values of t3 w /t q , the body perturbation becomes predominant 
compared with the heat perturbation. Therefore H c increases with M» 
for M ro > 1, a behavior directly opposite to that for a flat plate 
which does not offer any body perturbation (see fig. 8). On the other 
hand, for large values of ^ W /T 0 (see fig. 9 ), the body perturbation 
becomes less important compared with the heat perturbation. Therefore 
N c decreases with Moo for M* > 1, a behavior comparable with that for 

a flat plate (cf. fig. 6). 

From the above formula for the heat delivery, it is seen that Q 
is proportional to the coefficient X^P~ or X^Pr x Re. In the present 
theoretical investigation, this coefficient has been assumed as independ- 
ent of the temperature gradient in the neighborhood of the wire, -and 
the viscosity has been assumed negligible. In view of this simplification 
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and of the uncertainty attached to the value of X, the theoretical 
value of the coefficient, calculated for a constant temperature of 17° C, 
seems too high compared with the experimental ones. By comparing the 
formula of King with the experiments on cylindrical hot-wires at zero 
Mach number, it is seen that the ratio of experimental coefficient to 
theoretical coefficient is in the range of 1.5 to 1.7. These consider- 
ations have to be taken into account in the numerical interpretations of 
the results. 


Some Concluding Remarks on Heat Flow and 
Some Notes on Effects of Viscosity 

In this section various details of the reasoning are reviewed, and 
the effects of viscosity are studied. This will furnish at the same time 
an opportunity of pointing out briefly the various principles and hypoth- 
eses underlying the theory, and the difficulties which herewith are 
connected. 

In the present paper, it whs intended to elucidate some fundamental 
theoretical features of a heat flow, namely, 

(1) The generalized potential theory for a flow with rotation and 

variable heat energy 

( 2 ) The distribution of heat energy furnished by a heat source 

(3) The distribution of heat flux 

( 4) The amount of heat delivery 

(5) The perturbations of the flow by heat 

A simplified model was chosen for the study; it consisted of a 
straight array of continuous heat sources, placed parallel to a stream 
of uniform velocity. The heat sources on the array were of such an 
intensity distribution as to furnish a prescribed surface temperature, 
for instance, a uniform temperature. In order to fix the idea, such an 
array was called an "infinitely thin flat plate," which has the only 
significance of indicating the location of the heat sources giving a 
constant surface temperature, without offering any body perturbation to 
the flow, for example, in the form of a boundary layer. 

The theory of generalized potential flow has the advantage of per- 
mitting the use of the potential-streamline coordinates, and hence the 
application to more complicated shapes of heat source distribution, other 
than the above straight array. In fact the use of such coordinates can 
transform a complicated distribution into the straight array, which is 
therefore the most fundamental distribution. 
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The simplified model has served to illustrate certain physical 
features . However, it may, of course, deviate appreciably from an actual 
body supplying heat to the flow. For example, the flow in the boundary 
layer of the body should be actually retarded, and therefore the heat 
delivery should be less intense. But before going to such problems of 
a nonuniform stream, it was necessary first to study the more fundamental 
problem of a uni form stream, which gives a clearer analysis of the pre- 
dominant physical characteristics. The heat delivery by a blunt body, 
for example, a cylinder, gives rise to additional, difficulties which were 
not investigated, because the nonuniform flow behind a curved detached 
shock wave, as produced by the blunt body, is in itself also a difficult 
problem, even without heat introduction, and must be left out of the 
scope of the investigation at the present time. In that respect only 
some phenomena behind such a wave have been roughly stipulated in the 
preceding section, illustrating the difficulties which arise, rather than 
their explanations . 

The theory of the heat flow (nonadiabatic, rotational, and with 
variable total energy), in connection with the above enumerated aspects, 
was elucidated in the preceding sections by neglecting the effects of 
viscosity, a procedure generally adopted in the theories of such a com- 
plicated flow. In the following lines the effects of viscosity will be 
estimated. 

For a viscous flow, the equation of energy is (see reference 6) 


oc 53!. _ _ vo2m + G> 

P P Dt Dt “ 9 


where 


D_ 

Dt 


— + w • grad 
oti 


$ is the dissipation function (<X> Is used for potential perturbation, 
and the symbol $ is therefore chosen for the dissipation function) 
and all other symbols have been explained in the first section of the 
present paper . 

By substituting for p from the Efavier -Stokes equation 


P — = -grad p + r iv 2 ^ 
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where t] is the coefficient of viscosity. 


p 2_/ c T + 2 !) . |p = xyap + 

DtV P 2 St 


^T)WV^W + 


For a stationary flow. 


/ w 2 ^ 

pw • grad I CpT + —J 


= XV2t + ^wV^w + ® j (79) 


This partial differential equation deter min es the spatial distribution 
of the temperature T, or, following equation (l), the spatial distrib- 
ution of d(x,y), that is, the additional temperature furnished by the 
heat source. As was done in the section Small Perturbations, smal 1 
perturbations will be considered by keeping only terms in -3 of the 
first order, and hence it is possible to transform equation (79) into a 
linearized partial differential equation for 8. From the equation of 
continuity ( 5) , and from the relation of the density-temperature ratio, 
it is observed that cHrp/dxj is of the same order as and hence 

the second-order terms 

dw. cha. 

— 7^ (i = 1, 2, 3) 

OXj 0X2 


are negligible. 

The difference between equation (79) and the energy equation (6) 
used formerly lies in the presence of the last expression within paren- 
theses, which represents the effects of viscosity. 


Now, . 


w V 2 * = ^V^ 2 - (dwi/dxj) 2 



52 


MCA TN 2^36 


The last term of the right-hand member is a second-order term in d and 
is therefore negligible. So is $ too. Hence by neglecting second- 
order terms in the last expression between parentheses of equation (79), 


pw • grad ^CpT + = X^T + 


y2 

2 


= i 4„i + ^ — 
s l p ^ 2 


If the Erandtl number is written as Pr = rjCp jx = 1, 


w2\ c -nP 

V^fcpT + — j = w » grad |CpT + J; — 


w 2 


The density p can. be transformed into T, and the temperature T can 
be further transformed into d , using the formulas of transformations (2) 
and (3)- By leaving out second-order terms, after transformation there 
is obtained 


* ■ *' £ 


where 


1 _ 

1 "k-1 . TT -1 

g* = m* hU = gp ro 


(80) 


and g, h, Poc and U are defined as usual in the first section of 
the paper. Equation (80) is the partial differential equation for the 
spatial, distribution of d, by taking into account the effects of vis- 
cosity. It is seen that equation (80) is of the same form as the non- 
viscous distribution represented by equation (23), except the coeffi- 
cient 2g in equation (23) is replaced by 2g' in equation (80). 

The spatial distribution of the temperature d, the flux distrib- 
ution, and the amount of heat delivery follow from equation (80), exactly 
in the same forms as formulas (29), (37a), and (^3) with the new coeffi- 
cient 2g' as follows: 



MCA TN 2436 


53 


■d = e g,X K (ff'r)u 
2jt 4 r ^° 


q(x) = q(0)e 2g X (l - erf \j2g ' x) 


with q(0) = 2\d w 2g*/ii 00 ( according to formula (4l)) 


N = ^-({Qxg'z + l) 


(81) 


It is seen that the viscous distributions of $ and q are "rounded 
down, " that is , not so sharp as the nonviscous distributions, and that 
the heat delivery is diminished by the viscosity. The effects of the 
viscosity, which increases with the free -stream Mach n umb er, are thus 
to modify the magnitude of the physical quantities considered above, 
without changing the forme of their differential equations and their 
formulas . 

From formula (8l) for the heat delivery, it is seen that N is 
proportional to ^g 7 , and hence to \|u. Now the horizontal velocity U 
is subject to perturbations by the heat transfer. Those perturbations, 
of the order of •a, were assumed small, and as a first approximation they 
are neglected in the elucidation of formulas (43) and (8l) for the heat 
delivery in a nonviscous flow and a viscous flow, respectively. Later, 
as a second approximation, the flow perturbations by heat were calculated 
in the sections explaining flow perturbations by heat sources at subsonic 
and supersonic speed, by neglecting the viscosity. It is an easy matter 
to substitute the perturbed velocity in the formula of unperturbed heat 
delivery in order to find the perturbed heat delivery. Probably the 
viscosity will again have certain effects on the flow perturbations. 

But, as the flow perturbations play only a second-order part in the heat 
delivery, those effects will not be studied further here. 

Curves are plotted in figures for the case of a nonviscous flow, 
illustrating the general properties of a heat flow. The corresponding 
curves for a viscous flow are not added, in order not to make the figures 
too cumbersome. 

The effects of viscosity will reduce the slope of the curves for 
heat delivery in figures 6, 8, and 9, without changing the intercepts 
with the vertical axis. The slope of the curves in figure 6 will have 


i 
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to be multiplied by a factor that is, where C]_ is 

given by figures 3 and 4. The slope of the curves in figures 8 and 9 
will have to be multiplied by the factor 




where is given by figure 4 in terms of Mo,. 


National Bureau of Standards 

Washington, D. C., August 9 , 1950 
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Figure 4.- Mach number -functions at a supersonic speed. 




Figure 5.- Nonlinearity of heat delivery at a subsonic and a supersonic 
speed in the case of a flat plate. Broken lines represent corre- 
sponding linear relations obtained by neglecting heat perturbations. 







Figure 6.- Relations "between heat delivery, Peclet number, Mach number, 
and temperature difference in the case of a flat plate. The rela- 
tion for Moo = 1 and d v /T 0 / 0 is not plotted, as it is not 
covered by the elliptic and hyperbolic equations of perturbations. 
For lgL 3 rgQ Mach, nunbctrsj "fch.6 theory d.ocs not* pr6cLic*fc Tsori&tions of 
Nusselt number at low Peclet numbers. 
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Figure 7*- Nonlinearity of heat delivery at a subsonic an d a supersonic 
speed in the case of a cylinder. Broken lines represent correspondi ng 
linear relations obtained by neglecting heat perturbation. P c , Peclet 
number for a cylinder. 
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Figure 8.- Relation ‘between heat delivery, the Peclet number, Mach 
number, and temperature difference in the case of a cylinder for 
^/Tq = 0.05. The relation for M ra « 1 and fl w /T 0 / 0 ie not 

plotted, as it is not covered by the elliptic and hyperbolic equations 
of perturbations. 
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